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ABSTRACT 



Low speed / low manoeuvrability aircraft are currently quite susceptible to being killed 
in attacks by the ubiquitous infrared missiles. A theoretical analysis applied to an en- 
counter simulation seems to indicate that it is possible to use the infrared jammer to 
defeat second generation infrared missiles. The theoretical analysis of a simplified case 
of a conical scan reticle with frequency modulation jamming leads to expressions for the 
target's position, as seen by the missile seeker, under no-jamming and under infinitely- 
powerful-jamming conditions. The intermediate-power case is dealt with by numerical 
analysis for a selected, non-optimal situation, as the closed form solution is not imme- 
diately apparent. The analysis indicates successful jamming in the situation studied. In 
the scenario where the infrared missile is an almost continuous threat during the air- 
craft's flight, infrared jammers and low visual signature paints, and perhaps low infrared 
signature paints, are short-term solutions that are potentially useful in increasing the 
survivability of these aircraft by reducing their susceptibility to infrared missile kills. 
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time available, to ensure that the programs are free of computational and logic errors, 
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tional verification is at the risk of the user. 
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I. INTRODUCTION 



A. BACKGROUND 

When low speed f low manoeuvrability aircraft operate from bases which cannot 
be secured to a sufficiently great depth (more than several kilometres), they are vulner- 
able to surface-to-air missiles (SAM). As small, man-portable SAMs are quite readily 
available, they constitute a high probability threat. 

When out of the SAM engagement envelope, air-to-air missiles (AAM) remain as 
potential threats to these aircraft. 

B. SCOPE 

This thesis will address a particular aspect of this threat - the threat by infrared (IR) 
guided missiles. This is because infrared missiles are ubiquitous and rather successful, 
as can be seen in Table 1 on page 2 ((General Dynamics, 19SS), (Loral, 1989) and 
(Nicholas, 1988)) and Table 2 on page 3 (Loral, 1988, Slide 87329). l 

For this thesis, we will consider aircraft that are already deployed (as opposed to 
an aircraft still to be designed). Therefore the scope of survivability enhancement by 
vulnerability reduction is more limited than increasing survivability by susceptibility re- 
duction. 2 

Hence we will concentrate on the (hopefully) easier task of finding measures to reduce 
the susceptibility of these aircraft to being hit by an infrared missile. 

Due to U.S. Government restrictions on dissemination of security-sensitive infor- 
mation. this thesis is required to be unclassified. This constrains the thesis to a more 
generic treatment of the problem. 

C. GENERAL SCENARIO 

Tor the purposes of this thesis, the scenario is that the aircraft launches from a base 
that can only be secured to a depth of a few kilometers. i.c.,the aircraft is within the 
acquisition and launch envelope of threat missiles (especially easily-concealed man- 
portable SAMs) on takeoff as well as during the landing phase of flight. (An example 

1 Quoted from Moore D.,(OSD), 23rd IRIA IRIS 19S5 IRCM Symposium, p.l A-2, 1985. 

2 The terms used here are from Ball (1985). The probability of kill of an aircraft ( P k ) is the 
product of probability of hit by a threat weapon ( P,) and the probability of kill given that the air- 
craft is hit ( P k; ,). Susceptibility is associated with the probability of hit whereas vulnerability is 
probability of kill given a hit by the threat weapon. 
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Table 1. QUANTITIES OF MISSILES 



(Shows the extent of the threat.) 



Name 


Designation 


Origin 


Type 


Deployed 


Quantity 


Acrid 


AA -6 


USSR 


AAM 


7 


7 


Alamo 


A A- 10 


USSR 


A A M 


7 


7 


Anab 


A A- 3 


USSR 


AAM 


1972 


0(3) 


Apex 


A A -7 


USSR 


AAM 


1975 


•7 


Aphid 


AA -8 


USSR 


AAM 


1976 (?) 


7 


Archer 


AA-11 


USSR 


AAM 


7 


7 


Ash 


AA-5 


USSR 


AAM 


1961 


0(3) 


Atoll 


AA-2 


USSR 


AAM 


Early 1960s 


7 


Chaparral 


MIM-72 


USA 


SAM 


1979 


10 000 + 


Firestreak 


- 


UK 


AAM 


7 


7 


Gaskin 


SA-9 


USSR 


SAM 


196S 


7 


Gopher 


SA-13 


USSR 


SAM 


19S1 


200 


Grail 


SAG. SA-X-5 


USSR 


SAM 


1972 


50 000 + 


Gremlin 


SA-14 


USSR 


SAM 


7 


7 


1IN-5A C 


1 IX- 5 A C 


China 


SAM 


19SS 


7 


Kciko 


Type SI 


Japan 


SAM 


7 


7 


Kukri 


V3B 


S. Africa 


AAM 


7 


7 


Magic 


R.550 


France 


AAM 


19S5 


10 000 


Mica 


- 


France 


AAM 


7 


7 


Mistral 


- 


France 


AAM 


7 


7 


Piranha (?) 


MAA -1 


Brazil 


AAM 


19S9 


7 


Python 


- 


Israel 


AAM 


7 


7 


RB-72 


RB-72 


Sweden 


AAM 


7 


7 


Redeye 


FIM-43 


USA 


SAM 


1966 


7 


Redtop 


7 


UK 


AAM 


7 


7 


Shafrir 


7 


Israel 


AAM 


7 


7 


Sidewinder 


AIM-9L.M 


USA 


AAM 


1977.19S2 


150 000 


Stinger 


FIM-92A 


USA 


SAM 


1979 


15 000 + 


Stinger- POST 


7 


USA 


SAM 


19S2 


53 000 



Note: O(n) means of the order of magnitude of 10\ 
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Table 2. KILLS OF AIRCRAFT BY WEAPON TYPE 



(Note the large fraction of infrared missiles.) 



Type 


Kills 


Fraction 


Infrared missiles 


135+ to 146 + 


93% 


Radar missiles 


2 


1% 


Gunfire 


8 4- to 10 + 


6% 



of this aspect of the scenario is an Antonov An-26 or Antonov An-32 aircraft operating 
out of Kabul.) T he aircraft then climbs to its cruising J operating altitude, typically be- 
tween 5 and 30 thousand feet. While it is at medium altitude, the more likely infrared 
missile threat is that of the AAM. 

Additionally, it is assumed that the launch platform of the infrared missile can reach 
a position such that the aircraft is within the acquisition and or launch envelope of the 
missile. The aircraft, which is usually unarmed, is unassisted by other aircraft and lacks 
the speed to run away from an airborne threat. 

Since infrared missiles are fire-and-forget weapons, multiple infrared missiles can 
attack the same target. Therefore, this threat must also be handled. However, we will 
also assume that the multiple infrared missiles are of one type, during one attack. 
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II. THE AIRCRAFT CHARACTERISTICS 



A characterization of the aircraft in terms of the factors significant to the infrared 
missile countermeasure problem are as follows:- 

A. GENERAL INFRARED SIGNATURE FACTORS 

The infrared signature of an aircraft is dependent, to various degrees, on a number 
of factors. Various energy sources (and sinks) cause heating (and cooling) of the aircraft, 
and also reflect off aircraft surfaces. The energy is distributed over the aircraft primarily 
by conduction and radiation (and convection, in the case of the engine exhaust plume). 
The surface temperature and the surface spectral emissivitv together determine the 
emitted power. The incident energy, geometry and spectral reflectance determine the re- 
flected energy component of the infrared signature. A summary list of the various com- 
ponents follows ( 1 Iigby, 1 972). (Where available, typical values are also quoted 
(Loral. 19SS, Slide 91537).) 

1. Emission component energy sources/sinks 

• Engine (thermodynamic cycle) 

(Hot metal typically emits 1 k\V/sr/pm at >1.5 pm) 

(Engine exhaust shroud typically emits 100 W/sr/pm > 3 pm) 

• Engine exhaust plume (Typically emits 50 W/sr/pm at 3.8 pm to 4.8 pm) 

• Avionics/1 leat pumps. 

• Solar heating 

• Aerodynamic heating via adiabatic effects and friction viscous effects. Aerodyna- 
mic hcatine is not significant at the speed at which these aircraft fiv. ( < 400kts 
IAS). 

• Radiation exchange with atmosphere and via nocturnal cooling. 

• Tree stream conduction (Skin surface emission is typically 10 W/sr/pm at 8 pm to 
14 pm) 

2. Surface emissivitv 

Total emissivitv (as opposed to spectral emissivitv) is around 0.9 for the more normal 
painted surfaces. 

3. Reflection component energy sources 

In addition to the emission component energy sources (which are also reflection 
component energy sources), the other energy sources include:- 
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• Sun - Sunlight reflecting of! aircraft surfaces, particularly when specular reflection 
occurs, is known as glint. At typical values of 100 Wfsr/nm at < 3 n m, sunlight 
reflection is a significant contributor to the infrared signature, when it is present. 

• Sky - This is indirect sunlight, scattered by the atmospheric components (e.g, 
clouds). 

• Ground (earthshine) - Again this is indirect sunlight, i.e., sunlight reflected off the 
earth's surface or infrared radiation coming from the heated earth s surface. 

4. Reflectance 

These values are highly variable in the infrared bands for painted surfaces, such 
as these aircraft would have. 

5. Geometry/Structure 

While not directly a source of energy, the geometry of the aircraft does deter- 
mine how the infrared energy is spatially distributed. Geometrical effects include:- 

• Multiple reflections - The "energy packets" reflect off several surfaces before they 
reach the infrared detector. This redirects the infrared energy in complex ways. 

• Absorption/re-emission sites - Where a site is in relation to the energy source will 
determine whether it will absorb energy by radiation and/or convection. 

• Shielding - The energy source is masked by aircraft structures between the source 
and the infrared detector. 

• Shadowing - Part of a surface that appears darker because the energy source is 
shielded from that surface by intervening aircraft structures. 

• Engine exhaust plume shape and size - This is a function of the engine operating 
conditions and the airflow field around the aircraft. 

• Aerodynamic heating - This arises from adiabatic compresssion of air and from 
drag. Where this occurs, and the degree to which it occurs, is again determined by 
the shape, size and position of the aircraft structures. 

6. Atmospheric transmission 

The atmosphere is a selective absorber. (Cooper, 1988. p.1-6) In the infrared 

bands, there are two windows, one at 3 pint to 5 fxm and the other at 8 /tm to 14 /an 

waveband. Most of the detectors used in infrared missiles seekers arc designed to work 

in these windows. Most detectors are still usinc the 3 /urn to 5 n m waveband window. 

•c 

B. INFRARED SIGNATURE VS. PHASES OF FLIGHT/MISSION 

The takeoff and clinibout infrared signatures of an aircraft are significantly greater 
than during the cruise, on-station, descent or landing phases. The exhaust gas temper- 
ature varies by a factor of 2 or more between these phases of flight. 
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C. INFRARED SIGNATURE MODELS AND MEASUREMENTS 



Many infrared signature models exist. However, the more readily available ones 
only model the engine and only the jet-engine. e.g.,PIREPS (GEC, 1976) and ASDIR 
(ASD-WPAFB, 1975). There was no trace in the open literature of infrared models 
dealing with the turboprop engine, which typically power these aircraft. A few, (e.g, 
HIDE (Higby, 1972)) do consider the other factors such as skin surface emissions. Nei- 
ther infrared signature models nor measured signatures of aircraft were available for 
inclusion in this thesis. It should be noted that there is a lot of similarity of these models 
to the surface rendering models found in CAD, CAM applications. 

The true infrared signature of the particular aircraft, which we desire to protect, 
needs to be verified and or determined by actual measurement with appropriate 
radiometers. (In the analysis that follows, only the relative radiant intensity of the 
jammer with respect to the aircraft is required. However, the absolute power of the 
jammer can only be specified if the aircraft infrared signature values are known.) 

D. VULNERABILITY 

Since vulnerability of an aircraft is defined as the probability of kill, given a hit, the 
definition of "kill" strongly influences the value of P k . Aircraft kill would be when the 
aircraft is destroyed. However, if a kill is defined as denial of mission, then the P k may 
be larger, since electronic equipment, other support equipment and human operators are 
all susceptible to serious mission-hindering damage. 

E. MANOEUVRABILITY 

Low speed / low manoeuvrability aircraft are quite far, in flight performance, from 
an agile fighter aircraft. The climb rate of the low speed / low manoeuvrability aircraft 
is typically of the order of a few thousand feet a minute and their speed is in the low 
hundreds of knots. Acceleration normal to the aircraft's forward flight is limited to a few 
times gravity. This prohibits evasion of contact by acceleration and jinking. 
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III. THREATS AND THE ENVIRONMENT 



A. INFRARED MISSILE SYSTEM 

Here are some typical features of infrared missiles. 

1. Airframe and propulsion 

Usually the infrared missile has a cylindrical body with cruciform delta and/or 
trapezoidal wings and canards. The nose is blunt and transparent at the infrared 
waveband of interest. Propulsion for AAMs tends to be from single stage motors. SAMs, 
on the other hand, are often two-staged motors, with a booster to get the missile up to 
flying speed rapidly and then a sustainer motor for continuing the flight. The flight per- 
formance is characterized by the acceleration and velocity vectors attainable as a func- 
tion of time, and the maximum ranges this gives. 

2. Seeker 

a. Sensor 

Earlier generation detectors were made with lead sulphide. Indium 
antimonide is more widely used nowadays. Lead sulphide detectors work around the 3 
Mm waveband whereas the detectivity of cooled indium antimonide peal s at around 5 
Mm. This allows the cooled indium antimonide detectors to work with radiation from all 
aspects with respect to the aircraft. This is in contrast to lead sulphide which can only 
work with the hot engine parts and engine plume. (Cooper, 19SS,pp. 2-22, 7-21) 

To accomodate the wide range ( > 70 dB) in the scene and target irradiance 
as the missile flies towards its target, the detector is normally backed by an amplifier 
with Automatic Gain Control (AGC), so that the downstream processors have a nor- 
malized signal to work with. Significant parameters here are the attack and decay re- 
sponse times of the AGC and the saturation level. 

b. Optics and stabilization 

There is usually a telescope in front of the detector to gather the infrared 
energy. To maintain the seeker line-of-sight with narrow beamwidth optics, the telescope 
is normally decoupled from the missile body motion, i.e., the telescope is stabilized in 
inertial space. Significant parameters here are the maximum slew rates and the gimbal 
limits and the instantaneous field of view which is a function of the focal length of the 
optics and the size of the detector. 
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c. Target discrimination! position determination 

The determination of the target existence and position in the presence of 
background scene 'noise' is accomplished by various means. Early generation missiles 
had rotating reticles. In this scheme, information is extracted from the envelope of the 
detected signals. i.e v amplitude modulation (AM). This was followed by stationary 
reticles with the scene moved (nutated) over the reticle. Information, in this case, is 
carried in the frequency modulation (EM). Both these schemes transform spatial infor- 
mation into time and lose spatial discrimination. More recently, pseudo-imaging sensors 
have been used that scan the scene with a spot (or spots). For mechanical convenience, 
rosette scan patterns are used. The size of the spot is usually of the order of magnitude 
of the target image size. In the future, even more information extraction will be possible 
with imaging systems that can determine target features within the whole target. The 
effort here is to increase the ability of the sensor to determine smaller signature targets 
at longer ranges in the presence of higher scene 'noise', ( from both natural and inten- 
tional 'noise'.) 

3. Guidance 

The guidance law for many (perhaps most) infrared missiles nowadays is pro- 
portional navigation of some sort. The proportional navigation 'constant' may vary 
during the flight of the missile. 

The tracking loops are used to control the missile. The details of these deter- 
mine the missile response to aircraft manoeuvring and target position noise (e.g., glint). 
The usual tradeoff between controllability and stability applies. 

4. Steering 

Steering is accomplished by one or a combination of aerodynamic surfaces or 
vectored thrust. To save mechanical complexity, some missiles can only apply normal 
forces in one axis relative to their body coordinates, and depend on the rolling of the 
body to align the force in the direction desired. 

5. Fuze 

Most AAMs are both proximity fuzed and contact fuzed. The field of view of 
the fuze and the triggering range and velocity determine the detonation point. Radar and 
infrared (coherent and non-coherent) fuzes are employed. 

Smaller man-portable infrared SAMs are contact fuzed. Man-portable missiles 
need to be light in weight. Putting in a more complex, and hence larger, fuze, would re- 
quire a reduction in warhead size, for a given weight limit and missile maximum range 
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(amount of propellant). The smaller warhead would then be less effective in killing an 
aircraft at the longer range at which the fuze could trigger. 

6. Warhead 

Missile warheads are usually high explosive blast or high explosive 
blast fragmentation warheads. The pattern and the nature of the damage mechanism 
dictate how the aircraft will fare in a near miss. 

B. MISSILE ENGAGEMENT PHASES 

The engagement process goes through the phases of detection, acquisition, launch, 
flyout and detonation. Intervention in the engagement process is possible at any stage. 

C. LAUNCH PLATFORMS 

The launch platforms for SAMs range from a man to vehicles. They can be very 
difficult to locate. Usually AAMs are carried by fighter aircraft, though, of late, they 
are increasingly being carried by helicopters. 

D. OPERATIONS 

The method of operation of the missile, e.g., where the threat is likely to be, when, 
etc., are significant, but the combinations are many and they depend heavily on 
intelligence-derived information. This topic will not be discussed further here. 

E. ENVIRONMENT 

The urban environment that is given in the scenario means that the background 
clutter is higher than it would be in a rural environment. This makes missile engagements 
of targets more difficult in such environments. 
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IV. RESPONSE 



In the scenario of this thesis, the aircraft is highly susceptible to successful attack 
by infrared missiles. Countermeasures have therefore to be devised that can increase the 
likelihood of the aircraft continuing its mission. Some possible countermeasures are now 
discussed. 

A. LAUNCH SYSTEM DESTRUCTION 

Destruction of the infrared missile launcher is difficult because man-portable infra- 
red SAMs are easily concealed. 

B. LAUNCHER ACQUISITION ENVELOPE AVOIDANCE 

In the scenario given, avoidance of the acquisition envelope of the launcher is not 
possible due to the limited geographic depth assumed. 

C. LAUNCHER ACQUISITION ENVELOPE REDUCTION 

Reducing the acquisition envelope of the launcher is possible because the acquisi- 
tion is often visual for infrared missile engagements. Low visual signature paint schemes 
are required. 

D. THREAT PROPAGATOR DESTRUCTION 

Destruction of the missile in flight can be achieved either directly by directed energy 
weapons (DEW) or by triggering the fuze of proximity fuzed missiles. The DEW ap- 
proach involves first detecting the missile, tracking it, and then attacking the airframe 
or tracking/guidance system, or the fuzing system. DEWs (laser, electromagnetic or 
charged particle beam) are not yet available for operational use. Electromagnetic pulse 
triggering of the warhead is also a future option. Pointing of those weapons is also 
problematic because the missile has to be tracked first. If the spot scanning seeker proves 
difficult to avoid (see below) , then it may be necessary to make this a viable option. 
Proximity fuzes tend to be sensitive in the direction perpendicular to the longitudinal 
axis of the missile so that they are difficult to access for self-protection i.e.,from the 
forward direction. While it is possible to trigger them, it may not be straightforward to 
put enough signal into the fuze early enough to be useful. 
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E. THREAT PROPAGATOR AVOIDANCE 

The threat propagator can be avoided by out accelerating the missile and/or by 
manipulating the seeker's tracking ability. It is very unlikely that low speed / low 
manoeuvrability aircraft can out-manoeuvre a missile because these aircraft cannot 
generate a large enough acceleration vector and or displacement from the missile's flight 
path. Also, this method would require a reliable means of determining the presence and 
location of the missile. Manipulating the seeker can be achieved by reducing the signal 
provided to the seeker, increasing the noise at the seeker or providing other (false) tar- 
gets for the seeker. 

1. Reduce signal 

This can be accomplished by reducing the aircraft's infrared signature. This can 
be done intermittantly. e.g.,by engine power reductions at critical points in time, or 
continuously. For the latter, hot parts can be shielded (insulated) and airflow altered to 
mix with and or cool the hot parts. This will require structural and aerodynamic work 
to execute. 

For reflected energy reduction, low infrared signature paints can be considered. 
A potential problem with these paints is that the thermal load on the aircraft will in- 
crease (Ball. 19SS. p.293). 

Another technique is ■' v 'capturing' the Automatic Gain Control (AGC) of the 
amplifier of the missile's detector. Essentially, this involves depositing enough energy 
on the detector for a brief period (of the order of the attack time of the AGC) to cause 
the gain of the detector amplifiers to decrease to the extent that the real target signal is 
lost in the noise. The period of energy deposition has to be brief, or else the real target 
signature is actually enhanced, which is the opposite of the intention. To keep the de- 
tector blind for a prolonged period, this energy pulse is repeated before the AGC can 
decay completely. The overall intent is to keep the signal available to the missile proc- 
essing electronics as low as possible. The technique works particularly against missile 
seekers that derive information from the envelope of the signal, as the first generation 
AM seekers do. 

2. Increase noise 

Noise can be increased either by tactical use of natural sources, such as by 
placing the target up-sun of the missile launch platform, or by deliberate use of high- 
powered infrared sources. For the latter case, if the noise source does not saturate the 
detector, it may very well provide a much stronger signal for the missile to home on. i.e v 
become a beacon. 
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3. False signals 

Flares are a commonly used means of countering infrared missiles. They are, 
however, expendables. In the case of these aircraft which will remain (because they 
generally have low excess power) for several minutes in the acquisition envelope of the 
infrared missile, reliable warning of the presence of the missile must be available to 
trigger Hare release. Current false alarm rates and probability of detection of present 
day missile warning equipment are not adequate for this task. Because of the lack of 
manoeuvrability of these aircraft, separation of the flare from the aircraft, if needed, will 
not be aided greatly by manoeuvre. Rather, the ejection speed of the flare and gravity 
will be the determinant of the rate of shift of the centroid of the target and flare signa- 
tures.? Moreover, flares with the combination of the right spectral quality (to defeat 
multi-colour seekers), quick burn build up and long burn time with low weight, are not 
readily available. Flares may also be difficult to use on takeoff, if the takeoff path is over 
populated areas. 

Modulation of an infrared source to provide misleading signals to the seeker 
can be done. (Devices used for this purpose will be called infrared jammers.) Against 
AM seekers, in addition to the AGC capture, the phase of the envelope variations can 
be manipulated by the infrared jammer and so the position vector determined by the 
seeker will be false. We do not deal with these seekers here as it is assumed that they are 
relatively easy to deal with, 'flic second generation missile seekers tend to be conscan 
(FM) types. An analysis done to determine the potential of the infrared jammer to de- 
ceive such a seeker has been carried out and shows that it potentially can do so, as will 
be discussed later in Section “V. Conscan (FM) seeker vs. infrared jammer" of this 
thesis. Because of the narrow instantaneous field of view of a spot scanning seeker, it is 
unlikely that significant false signals can be inserted into these seekers in this way. 4 



3 It may, on occasion, be useful to have a low separation speed against a spot seeker missile 
that uses kinematic lolloping to try to defeat flares. In that case, the missile miss distance will be 
small. Therefore, this can only work if' the missile does not have a proximity fuze. 

4 While the spot seeker missile must be dealt with, readily available information on these 
missiles is lacking. Therefore, with the time constraints of this study, we can only speculate on 
possible countermeasures against these missiles. 
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V. CONSCAN (FM) SEEKER VS. INFRARED JAMMER 



A. BACKGROUND TO THE ANALYSIS 

This is a reticle in which the image is moved (nutated) over a stationary reticle. The 
resulting frequency modulation of the signal out of the detector can be processed to ex- 
tract the 2-D position of the target relative to the reticle centre. 

The following analysis is derived from papers by Gedance( 1961) and Suzuki( 1979). 

Gedance's analysis is for a target in the form of an illuminated disk of uniform and 
constant irradiance, and a spoked (or wagon wheel) reticle with a transmissivity function 
that is either one or zero as a function of the azimuth of the reticle. His analysis stops 
at the determination of the instantaneous frequency generated by the target image on 
the reticle. lie does not derive the 2D position of the target. 

Suzuki considers a point target for a reticle similar to that analysed by Gcdancc. 
Suzuki then shows the relations necessary to obtain the 2D position of the target. 

Neither Gedance nor Suzuki deals with the jamming case which the following 
analysis does include. 

B. ANALYSIS 

In the following analysis, we want to determine the apparent target position, as seen 
by the seeker, while under jamming. To do so, we obtain, in order, the flux seen by the 
detector, then the instantaneous frequency of the flux signal and, finally, the expressions 
for the target position, as determined by the seeker. The expressions for the target posi- 
tion will be for the two extreme cases of no-jammer and for an infinitely-powerful- 
jammer. Only one case is dealt with here for the intermediatc-jammer-power range. This 
is because the general solution to this case is not immediately apparent and hence a 
numerical analysis had to be used. 

C. FLUX SIGNAL 

Assume a uniformly illuminated target disk of radius 5 moving in a circular path 
of radius r about the origin 0 ROT (See Figure 1 on page 14.) The position of this origin, 
with respect to the reticle centre 0 RET , is in fact the "position" of the target with respect 
to the seeker boresight. The target image centre is at ( D , a) with respect to the reticle 
centre 0^^ j- . 
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Target Image 




Orot 



Figure 1. Reticle-Target Geometry 
In general, the flux, F(t) , transmitted through a reticle is given by 

/ (0 = /co f T(s,e)dA ( 1 ) 

J A 

where 

I(t) = irradiance of the image, 

A = area of the image, 

(5, 0) = polar coordinates relative to the reticle centre 0 RET , 

T(s, 6) = transmissivity of the reticle. 

Ihe equation of the circle that describes the target image is given by 

S 7 = s 2 + D 2 — 2 sD cos(a — 0) (2) 
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Solving for 5 , in the above equation, we get the two roots 



s x = D cos(a — 9) — y' D 2 cos 2 (a — 9) — (D“ — <5‘) (3) 

s 2 = D cos(a - 9) + J D 2 cos 2 (a - Q) - ( D 2 - <5 2 ) (4) 



The elemental area, dA , is approximated by a rectangle of width 
side lengths (s 2 — 5,) . Therefore, the incremental area is 




and 



dA 




A + -h 
2 




( 5 ) 



From Eqns. (1), (3), (4) and (5), the flux F(t) is now 



’2+ a res in 



no - no 20 



Hi). 

) 



Tls, D) coija - ()) ' I.r cos'(a - 6) - (D‘ - 6‘) d6 (6) 



if the reticle transmittance is independent of the radius, s , i.e v the reticle is a radially 
symmetric reticle, then 

<■. 

T(s, 9) = T{6) (7) 

And if T{9) is a periodic pattern of period 2 L ( = ) , it can be described as a Fourier 

Series thus: 

GO 

T{9) = ~Y+ cos^ ■ nJ Y~ + <f>^j 

n=] ( 8 ) 

oo 

c 0 

= — + cos (nm9 + 4>„) 

n= I 

where 
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n me 



-jnmO 



dO 



( 9 ) 



Inserting this description of the reticle transmittance function (Eqns. (7), (8) and (9)) in 
the flux expression (Eqn. (6)), and reversing the order of summation and integration, 
we can write 



m = /(/) 



Aq + 






(10) 



where 



A 0 = 2D 



2 + arcs in 



(i)c 0 



-j- cos(a - 0\/D 2 cos 2 (a - 0) - ( D i - <5 Z ) dO (1 1) 



2 



<i) 



and 



• k + arcstn 



A 



Ki) 



/, A 2 _ 1, 



„(t) = 2d\ K D J c n cos (nmO + <f> n ) cos(* - 0 )Jd 2 cos 2 (a - 0) - {D 2 - 6 2 ) dd (12) 

J . ( 6 \ 

3f— — j 



Evaluating A 0 (Eqn. (1 1 )). we get 



A 0 — DCq 



* *+ arc sin 



•a) 

a) 



cos (a - 0\/d 2 cos 2 (a - 0) - (Z) 2 - <5 2 ) </0 (13) 



Let O' = a — 6 . Then get 



Aq — Dc 0 



csin (-0 ) 

<i) 



cos(0')^/Z) 2 cosV - (Z> 2 - (5 2 ) 4 - O') 



(14) 
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Since the integrand is even, 



A 0 = 2DSc 0 



= 2Ddc n 



S' a res in 
f*arcs in 



ln {-o) 



cos (d')J[ -J-) 2 C0S 2 $' - { -yV+l dd' 



(i) 



(15) 



cos(0'y 1 “ (-y) 2 sin20 ' 



From tables (Beyer, 1987), Eqn. (15) becomes 
A 0 = 2D6q 



sin 



O' I 



V 1 - ( -f ) 2 sulV + 1 77 arcji ”( 7T sin °') 



arcsin 
0 



KD 



(16) 



0 A 7 
= — nd 



Define 



. A i c 0 -2 

— Aq — 



(17) 



Likewise, for A„ (Eqn. (12)), 

i+ arcsin ^ ^ 

arcsin ^ -j- ) 



A n (0 = 2Dc,J 



cos(nm0 + 4 > n ) cos(a — 0)-Jd 2 cos 2 (« — 6) — (D i — <5 2 ) dO (18) 



For simplification, change the integration variable to O' = a — 0 . Also, expand the first 
cosine term in Eqn. (IS) . 'I hen 



A„(0 = 2 Dc n 



= 2 DSc„ 



— arcsin ( -j- ) 

rr " iM (-zr) 
rcsin ( ) 

ircsin^jj) 



cos {nm(ct — 9') + <f> n ) cos O' 6 /( — ) 2 cos “0' — (( -^r- ) 2 + 1) </( — 0') 

Vo o 



cos(mm + (/>„) cos cos 0'^yl - ^ -y- ^ 2 sin 2 #' 



(19) 



+ sin(/a?!:< + <f> n ) sin nmO' cos 0'^ 1 — ( ■— j 2 sin 2 0' 



dO' 
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The first term in the integrand of Eqn. (19) is an even function, while the second term 
is an odd function in O' and the integral of an odd function taken symmetrically about 
the origin vanishes. Therefore, 



f« 

A n (t) = 2 Doc n cos (nrnx + <£„) v " '2 cos 6' cos nmO' j 1 

.It can be shown (Gedance. 1961, Eqn. 9 and 10)5 that if -j- <|1 and ~jj - < 1 , then 



Kl). 



^ - J- V sin 2 0' dir (20) 



nm 






2 cos 0' cos nmO'J 1 - ^ -j- ^ 2 sinV '*7, ^ J11 ~ ^ (21) 



where J, 



ijrnd 

D 



is the first order Bessel function of the argument 



( nmS \ 

D )■ 



Then 



A w (r)as ft m D3c n cos{mntx + 4> n )J l 



timd 

D 



( 22 ) 



From the geometry of the situation (See Figure 1 on page 14.), the triangle formed by 
the centres of the reticle (0 RET ). target and rotation of the target image about the reticle 

ROT ) C $ 

D 7 = R 2 T + r 1 — 2Rjr cos(o) r r — 0 T ) (23) 



and 



r sin cop — R r sin 0 7 

a = arcian\ — 

r cos ay — R T cos 0-j 

Establishing p as a normalized offset distance, — ■ , then 



I) = r^J 1 — 2 p cos(c at —0 T ) + p 



(24) 



(25) 



and 



a 



arctam 



sin a ) r t — p sin 6 r 
cos cop — p cos 0 T 



( 26 ) 



5 See Appendix A for a demonstration, by series expansion, of this approximation. 
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A„ can then be written in the condensed form 



A„(r)*A„ cos 0„(O 



(27) 



where 



A„ = 



r 



nm V 



rJ 1 — 2 p cos(cot — 0 T ) + p dc n JA 



nm5 



■-A 



2 p cos(a»/ — Q-j) + p 



(28) 



and 



A / sin ay — p sin 0 T \ 

Q n {r) = nm arctanl — + <f> 

\ cos (Jj r t- P cos Q r J n 



(29) 



Note that for , A n (r constant. Note also from Eqns. (27) and (29) that the infor- 
mation on the position (p , d T ) is carried in the angle, i.e. the target image motion and 
reticle interaction gives rise to a frequency modulation (FM). If the reticle pattern is 
such that the conditions —■ <^1 and — < 1 can be met6 , then 

F(0*/(t)(^o + Y,A n {t) cos 0„(O^ (30) 

where A' is large but not violating the condition. 

D. INSTANTANEOUS FREQUENCY 

We now develop an expression for the instantaneous frequency of the flux signal 
as a step towards determining the information carried in the EM. 

The flux signal. F(i) is then bandpass filtered. Let's call the bandpass filtered signal 
F b {t). In general, a bandpass filtered signal can be represented in quadrature form or 
complex envelope form, i.e., 

F b (t) = P(t) cos c o c t - <2(0 sin a> c t (31) 



6 A seeker with a wagon wheel reticle would normally be designed to have the target image 
of approximately the same size as the width of the transparent spokes of the reticle, at the radius 
at which the target image rotates on the reticle. This would give the greatest targct-to-background 
discrimination. Also, during a subsequent bandpass filtering operation (which will be describe later), 
the higher order components (associated with n in one of the conditions) of the reticle-generated 
signal are removed. Thus, in practice, in the vicinity of the design point of the reticle seeker, the two 

conditions of ~ 1 and < 1 can usually be met. 
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or 



F b {i)^Re[g{t)^ (32) 

where, in this instance, if in is the number of cycles of transmittance changes that a point 
on the image goes through in one rotation about the reticle, then oj c = moj, . The two 
forms given in Eqns. (31) and (32) are related by 

g(‘) = P(')+jQ(‘) ' (33) 

From the quadrature components P(i) and Q(i), we get the instantaneous frequency 
of that signal' (Roberts, 1977) as 



COi(l) = co c + 



P(0Q(0-Q(m0 

r\,) + q\,) 



(34) 



In a step reminiscent of synchronous or coherent demodulation of an FM signal, but 
using the scene rotation frequency as the reference frequency, we can extract the' -in- 
phase and quadrature-phase components of the signal. Later, we will show that this in 
fact gives the position of the target. To get an easily recognizable expression for the 
position in the end, we also include a scaling factor ( ) , at this stage. The in-phase 
component, /( r) , and quadrature phase component, q( t) , are described by 












COS C O r l 



sin u) r i j 



di 



(35) 



Note that we can only find the position vector a half-cycle behind current time because 
the integration takes place over a full cycle. 

E. SIMPLIFIED RETICLE AND JAMMING SIGNAL 

It is rather tedious to continue the analysis generally. So that we do not lose sight 
of the forest for the trees, we shall consider a special case of a reticle and of a jamming 
signal. 



7 Essentially the instantaneous frequency is given by the differentiation of the phase of F b {l) , 
with respect to time, and offset by the carrier frequency. 
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We can simplify the reticle transmission function to a single sinusoidal azimuthal 
variation, reducing the Fourier summation description of the reticle to a single term, 
besides the constant term. i.e v 



T{6) = y + 4" cos[m(0 - <£,.)] 



(36) 



In terms of the quantities defined in Eqn. (8), 



c 0 = 1 



c \- -) 



c n = 0 for n > 1 

Then, from Eqns. (17), (28) and (29), we get 



An ~ 



710 



(37) 



A](t) = •— Or N / 1 — 2 p cos(c oi —0 T ) + p z 



mS 



l 2 

Ty- 1 — 2 p cos(aj/ k —d T ) + p 



(38) 



and 



0,(0 = m 



arctanl 



sin oj r t — p sin Q T 
cos ix> r t — p cos d T 



~ <i> r 



(39) 



We now have a description of the component of the signal that results from the inter- 
action of the target and the reticle, in the seeker. 

F. VARIABLE IRRADIANCE FROM JAMMING 

To investigate the potential for jamming, we apply an optical signal combining a con- 
stant target irradiance with a sinusoidal jamming signal. The irradiance , I(t) , is then 
described by 



AO — Itgt + (/C 1 + cos 0/0] (40) 

where I, g , is the target irradiance 

and /[I + cos 0 ; (O3 is the jammer irradiance which oscillates between 0 and 2 1 ; . 
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Rearranging Eqn. (40) , we get 



7(0 = I c + Jj cos 0/r) 



where 



is the unvarying component of the irradiance. 

The flux signal for this special case is 

T(0 = 1 I c + Ij cos ©/OIOo + cos 0i(OD 

After high-pass filtering, the D.C. term, I C A 0 , will be rejected by the filter, 
bandpass signal 



4 

W = y p k cos <t>k(r) 

k=\ 



where 

F i = IM 0 0, = 0,(0 

r=i/i o <t>2 = 0.(0 

fj = W) 03 = © y (0 -0,(0 

^= 7^(0 04 = 0/0 + 0.(0 

In complex envelope form, 

F/0 = 7?c^(/)^} 



where 



? (0 = Vf/' , ' ) - w ' ,) 

>!=1 



For compactness in subsequent expressions, let 

<J>*(O = 0*(O-«V • 



(41) 

(42) 

(43) 

and we get the 

(44) 



(45) 

(46) 

(47) 
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To find the instantaneous frequency of the bandpass signal, we need first to find the 
quadrature components ( P{i ) and Q(t)) of the bandpass signal. This is done using Eqns. 
(31) and (32) with Eqns. (46) and (47). We then get 

fv> = kv, sm) 

-In «.«*) ,4S) 

n=\ 



and 



Also. 



( 2(0 = 0 } 

4 

= sin <D a (0 

n - 1 



(49) 



4 4 

t 2 (0 + Q\i) = Z Z ^ C0S (^W - <*>/«) (50) 

k= 1 /=1 

Taking the derivatives of the quadrature components described by Eqns. (48) and (49) , 
we get 



P{t) = 



sin <D*(0 



n = 1 



(51) 



and 



4 

(2(0 = ^^(O* 7 * cos <D*(0 
«=1 



(52) 



Also, combining Eqns. (48), (49), (51), (52), gives 



4 4 



P(t)Q(t) - Q(‘)H‘) = Z Z <*>*(0 V, cos(O,(0 - 0,(0) (53) 



k= 1 /=! 
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Substituting Eqns. (53) and (50) in Eqn. (34) , the instantaneous frequency in 
is described by 



4 4 

XZ <S> k {t)F k F t cos(c]>,(/)- 0,(/)) 

wXO = + — r 4 

£ cos(<D*(0 - <DXO) 

Ar=1 /=1 



Now we define a jam-to-signal-plus-jam ratio, /? y , and further restrict the analy 
case of small target offset from centre. Let Pj = ~r = ~r~ ~ y • If P< 1 ( 

*C *Tgt * V 

target offset from centre), then 



A 0 *Aft)*A c 



Hence 



h^c Fj Ja_ Ja_ n 
1 C A C * F , * F, * /•) 

and we then get 



EE^W/CosWO-^O) 

“KO = w c + 

JE^/COS^W-OKO) 



where 



and 



f I if k = 1 
^ = // A * 1 




«/ /=! 
if l* l 



this case 



(54) 

sis to the 
i.e v small 

(55) 



(56) 



(57) 



(58) 



( 59 ) 



24 



G. SPECIAL CASES 

Let us now examine this result in the two limits of /?, = 0, no jamming, and /? ; = 1, 
the infinite jamming case. 

1. No-jamming case 

Consider the no-jamming case when f = 0 . i.e v when /?, = 0 . Then 




if k = 1 
if k* 1 



( 60 ) 



and 



P,= 




if /= 1 

// /#! 



(61) 



Substituting Eqns. (60) and (61) in Eqn. (57), we get 



+ 



(!),(/)( If 



or 

- + </>,(>) - 

= G,(r) 

1 — p cos(a v — 67 -) 



(62) 



= maU 



1 — 2p cos {(o r t — Of) + p 



Now that we have the instantaneous frequency of the signal, we will show that the co- 
herent demodulation of Eqn. (35) does indeed give the position of the target, f'or the 
in-phase component. 



I 

nm 



*r+« 

CO, 

c Oj([)r cos c o r idt 

J 0 j — K 



CO, 




1 — p cos(o) r / — Of 
1 — 2 p cos(a ) r t — Of) + p 1 



cos oj r tdt 



(63) 



“V 

Defining the change of variable O' — ay — 0, and expanding the trigonometric terms, we 
convert the integral into the form, 
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( 64 ) 









cos O' cos 0 T — sin O' sin 0 T — p cos 0 T cos O' + p sin #7- sin O' cos O' 



2 r cos 



1 — 2p cos(a> r f — 0 T ) + p~ 
cos O' — (p'2) cos 2 O' + (p '2) 



rf0' 



1 — 2p co$(co r t — 0 j) + p 



■dO' 



From the tables of integrals (Gradshteyn 1980) we can evaluate this definite integral, as 






2r cos 0 ^ 



[f] 



= rp cos 0 T 



(65) 



— RfCos Of 



which is in fact the .v component of the target position. Similarly, from the 
quadrature-phase component, we get 



cj — Rf sin 0 f (^^) 

which is the y component of the target position. 

Hence the target location in two dimensions is now determined. (Note that 
whether or not there is jamming, the missile seeker, not knowing any better, will still 
perform the same procedure, described above, to try to extract the target position.) 

2. Infinitely powerful jammer 

Let us now consider the case in which we have an infinitely powerful jammer, 
i.e. /?, = 1 This case gives an instantaneous frequency expression that is not simple. To 
gain further insight, we will restrict the bandpass filter to an even narrower band about 
the carrier frequency ( co c ) of the flux-derived signal. The upper limit of the bandpass is 
just above twice the carrier frequency. Mathematically, we make d> 3 and <J> 4 zero. Then 
I2qn. (57) becomes 



co,(/) = co c + 



0,(r) + <i> 2 (/) + [<i>,(r) + ci> 2 (/)] cos[(I>,(/) - ch 2 (r)] 
2 + 2 cosC(I>,(r) — d> 2 (/)] 



- q- [0,w + 0/o] 



(67) 
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And if the jamming waveform is similar in form to the reticle modulation waveform (as 
given in Eqn. (39)), and the jammer's equivalent of the reticle rotation frequency is 
co rJ — co r , then the position coordinates of the apparent target location are : 



— (R r cos 8 r + Rj cos Oj) 
\ {R T s\n 0 T + Rj sin Oj) 



( 68 ) 



where Rj & Oj are the jamming parameters corresponding to R T & 0 T of the reticle. Note 
the effective pointing direction of the seeker to the apparent target is that for the 
centroid of the combination of the target signal and the jamming signal. 

H. SOME OBSERVATIONS ON THE ANALYSIS RESULTS 

The jamming signal cannot dominate the target signal because the jamming signal 
is itself modulated by the relative motion of the jamming spot and the reticle. 

The analysis here is of a reticle that has a transmissivity function that varies grad- 
ually in a sinusoidal manner in the azimuth direction. For manufacturing ease, most of 
these reticles in real missiles will have more sharply changing transmissivities. However, 
to limit noise, the signal from the detector will eventually have to be bandlimited. Thus 
this analysis is a good first order approximation for the real missile. 

It should be noted from the analysis that the conscan reticle seeker does not derive 
its information from the amplitude envelope. This being the case, AGC jamming tech- 
niques may not prove to be particularly effective against the conscan reticle seeker. 

I. SIMULATION 

Numerical analysis was performed to examine the intermediate cases of /?, . For 
simplicity and computation speed, only the narrow band version of the analysis was in- 
corporated in this simulation program, i.e.. 






2 2 



*>/(0 = w c + 



2 2 



Z Z M cos(<D*(0 - <t>fa)) 



( 69 ) 



= ma\ + 



0,(0 + /?>,(/) + [0,(0 + <j) 2 (f)]/?y cos[O,(0 - Q 2 (Q] 
1 + ^ + 2/? y .cos[O ] (0-O 2 (/): 
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Eqn. (69) was then inserted into Eqn. (63) and the corresponding equation for the y-axis, 
and was transformed into computer code. At the same time, to incorporate the missile 
aerodynamics, the 2D position vector computation from the seeker was inserted into a 
version of the TACTICS IV program (Meisberger,19S3) which was ported to a micro- 
computer. TACTICS IV is a time step simulation of a missile flight. It computes the 
relative positions and velocities of a target and a missile, takes the seeker-derived posi- 
tion (with noise added, if desired), applies a proportional navigation guidance law, and 
integrates the effects of a commanded acceleration to obtain the new position and ve- 
locity of the missile. The modified TACTICS IV, which we call TACTICS IV (IR) in- 
corporates the infrared seeker behaviour with an infrared jammer. The seeker routine 
was patched so that it will branch to the new infrared seeker routine when desired. 8 The 
parameters required for the seeker are essentially the parameters of Eqn. (69). In addi- 
tion, the tabulated output file of the missile and target flight profiles was modified 
slightly to be more easily read by a graphics program.9 The listing of the infrared seeker 
portion of that program can be found in Appendix B. 

An example of the result of an air-to-air missile flight simulation is in Figure 2 on 
page 29 . In this example, the missile is launched at the target from a position directly 
behind the target. The target does not evade, but continues in a straight and level flight 
path. Figure 3 on page 30 is the input data to the TACTICS IV (IR) program, for this 
example. In this particular case, the missile departs from the target in a spiral-like path. 
Note the increase in miss distance to 700 feet with the infrared jammer on, as compared 
to the 3 feet case when there is no jammer. 10 

This study should be continued with sensitivity analysis studies and optimizations 
to identify the optimal jamming parameters. 



8 Tills is done by making the 1SEEK parameter of the input data table equal to 5. 

9 The postprocesser described by Meisberger (1983) was not ported to the microcomputer. 
Instead, a commercial package. STATGRAPHICS, was used for the graphical output. 

10 The no-jammer case was simulated in a separate run, with /y effectively zero. 
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Missile Flight Profile 




Figure 2. 



Missile flight profile (\>ith jamming): 
left to right. 
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Missile and target flying from 
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Figure 3. TACTICS IV (IR) Input data for missile flight (with jamming) 
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VI. RECOMMENDED COUNTERMEASURES 



In summary, given the scenario considered here, the most promising short term 
countermeasure solution is the infrared jammer. The infrared jammer has a good chance 
of working if it has the following characteristics:- 

• For all aspect coverage, the infrared energy sources should radiate all around the 
aircraft. This would imply more than one source element, as there is no single point 
on the aircraft that has equal visibility in all directions. 

• The power output of the jammer source, in the band that the missile detector is 
using, needs to be at least equal to the target signature in that band and in the same 
direction. 

• The jammer should be able to modulate its infrared power emitted such that the 
waveform is similar to what the seeker forms when it views a steadily radiating 
target. 

The analysis done here is only a paper study, using what are intended to be (hope- 
fully) typical parameter values. The actual signatures of these aircraft under various 
conditions must be measured to provide accurate input. The infrared jammer source and 
modulation scheme must actually be verified by test and evaluation with real hardware 
to determine if the assumptions and analysis arc correct. e.g.,can we find IR sources that 
can provide enough in-band energy on the seeker's detector. We have not determined 
yet the optimal modulation scheme against the conscan (FM) seeker. Neither have we 
considered the optimal modulation scheme when there are more than one type of missile. 
Evaluation must be extended to these cases to ensure that modulation schemes, while 
effective against a particular type of missile seeker, do not decrease the miss-distance for 
other seeker types. 

This study must also be extended to the third generation scanning spot seeker 
threat, which is increasing with time. 
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VII. CONCLUSION 



Low speed / low manoeuvrability aircraft are currently quite susceptible to being killed 
in attacks by the ubiquitous infrared missiles. A theoretical analysis applied to an en- 
counter simulation seems to indicate that it is possible to use the infrared jammer to 
defeat second generation infrared missiles. The theoretical analysis of a simplified case 
of a conical scan reticle with frequency modulation jamming leads to expressions for the 
target's position, as seen by the missile seeker, under no-jamming and under infinitely- 
powerful-jamming conditions. The intermediate-power case is dealt with by numerical 
analysis for a selected, non-optimal situation, as the closed form solution is not imme- 
diately apparent. The analysis indicates successful jamming in the situation studied. In 
the scenario where the infrared missile is an almost continuous threat during the air- 
craft's flight, infrared jammers and low visual signature paints, and perhaps low infrared 
signature paints, are short-term solutions that are potentially useful in increasing the 
survivability of these aircraft by reducing their susceptibility to infrared missile kills. 
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APPENDIX A. VALIDATION OF THE APPROXIMATION 



S YIW\ 5 

We want to demonstrate that if -jy <^1 and ■ ^ < 1 , then 

run f. / D V • Z M , ,( 

— — J 2 cos 6 cos nmd 1 — ^ — j sin 6 dd 



nmd 

D 



( 70 ) 



For convenience, let 



A ( nmd 

X - ' D 



(71) 



and call the left hand side of Eqn. (70) ,f cm { x ) ■ Now the first order Bessel function in 
series form is 



J iM = 



_Y 
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oo 



k = 0 



/ A\k 7k 
( - 1 ) * 

2 2k k\(k+ 1)! 



(72) 



Expanding the first few terms of this series, we get 

aw* "I - * - iV * 3 + IsT * s ~- mi 

Expanding the integrand in the left hand side of Eqn. (70), using the Taylor series, and 
integrating, we get 



/cos(-V): 
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(74) 



If J<^1 , the terms in the second row of Eqn. (74) are small in relation to those in 
the first row and hence can be neglected. Computing the coefficients of a , 
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(75) 
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■Comparing Eqn. (73) and Eqn. (75), if 



nmS 

D 



< 1 , then 



/cosM^lM 



(76) 



which is what we set out to show. 
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APPENDIX B. TACTICS IV (IR) PARTIAL LISTING 





Program listing of the infrared seeker 
portion of TACTICS IV (IR). 

Requires TurboPascal 5.0 and the 
TurboPascal Numerical Toolbox. 

* ********..* ) 

unit tac8s; 

interface 
uses tac 89 ! 

{$1 f loatdef . inc> 

procedure cs j_seeker (time : float; 

var seeker_theta_prev . 
seeker_phi_prev . 
seeker_theta_cur . 
seeke r_phi_cur : float);* 

procedure GetOutputFi le_Name ( var OutFile : text* 

var FileName ; string; 
default_f i lename : string); 

procedure get_root_name_of_f iles(var FileName : string; 

de f ault_f ilename : string); 



implementation 



< $ I - } { Disable I/O error trapping } 

(SR*) { Enable range checking ) 

uses 

Integrat* Dos. Crt. Common. Funct; 

const 
m = 12; 
m_J = 12; 
r = 1.0; 
u_r_r e V = AO: 

dellw_r_rev_def ault = 0.39; 

phi c_r = 0.0; 

phi c_3 = 0.785A0: 

rho_Tnit = 0.0; 

rho_3_def ault = 0.A99; 

the t ac_T inlt = 0.0; 

the tac_ j default = 0.0; 

del_theta_ default = 3 . 1A 159265/3 ; 

beta_3_def ault = 1. 0/2.0; 

f ilt_cons = 1.0; 

max_move_step = 0.1; 

f ilename_cs j_de*ault = ‘CSJ.OUT'; 

focal_lengtn = 57.0; 

var 

w_r , w_r_j. del_w_c. del_w_r. beta_j ; float 
*_tj. y_t3* range_t 9 t. theta_tgt : float: 
x_t 1 . y_tl. range_tgtl- theta_tgtl ; float; 
var t eye. t_cyc 3 ; float; 
refit. dellt_2 - : float; 
u_l 1 m . l_lim : float; 
del_theta ; float; 
rho_j. rho_3_start ; float; 
del_w_r_rev T float; 

i ; integer; 

consl : text: 
rho- rho_step ; float; 
thetac_T7 thetac_3 : float; 
temp_beta : float; 
new_x. new_y : float; 
del_range. del_x. del_y ; float; 

F i 1 eName : string; 



LowerLimit. UoperLimit ; float; { 
Tolerance : float ; { 

Kaxlntervals : Integer: { 

{ to approximate the integral } 
Integral ; float; ( 

Numlntervals : integer; { 

{ to approximate integral } 

Error ; byte; ( 



Limits of integration ) 

Tolerance in the answer > 

Maximum number of subintervals used > 

Value of the integral ) 

Number of subintervals used } 

Flags if something went wrong } 



<$F + ) 

function beta(t i float) : float; 

var tt ; float : 

begin 

beta := beta_j; 

end: (functlonbeta) 



function theta_o(t ; float) : float; 

var numer, denom : float; 

begin 

numer :* s 1 n ( w_r« t )- ( r ho*s i n ( t he t ac_T ) ) ; 
denom := c os ( w^r " t )- ( r ho * c os ( t he t ac_T ) ) : 
tbeta_o ;s m** ( arctan ( numer / d e nom )-phl c_ r ) ; 
end ; 

function theta_j(t : float) ; float; 
begin 
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theta_ j 
end ; 



m_j*(ercten((sln( (w_r *de l_w r)"t)-Crho_j>'sln(thet0C..J)))/ 

(cos( (u_r*del”w~r )*t)-Crho_j"cos(thetac_j))))-phlc_j) 



function thetadot_oCt t float) : float: 
begin 

the tadot_o := C C l-rho«cos(w_r»<t-thet 0 C_T ) )/ 

( l-2*»rho»cos(w_r*t-thetac_T)*rho»'rho ) ) ; 

end ; 



function thetadot_j(t i float) : float; 
var numer, denom « float; 
be 9l n 

numer := C l-rho_j "cos ( (w_r*del_w_r ) * t-thetac_j ) ) ; 
denom : = ( l-2* , rho_3 “C osC (u_r» t*del_w_r )-thetac_ j ) + rho_j «rho_ j ) ; 
thetadot_J : = m_J *• ( w_r + del_w_r )*■ ( numer/denom ) ; 
end ; 



function phl_ j ( t ; float) » float; 
var tj * floit: 
be 9 l n 

tj := the ta_ j ( t ) ; 
phl_J t j-(m*w_r*t ) ; 
end ; 

function phldot_j(t « float) : float; 

var tdj : float! 

begin 

tdj := the tadot_j ( t ) ; 
phldot_j j= td j- (m«w_r ) ; 
end ; 

function phl_o(t : float) : float; 

var tho : float: 

begin 

tho := theta_o(t); 
ohl_o i = tho- (m*w_r * t ) ; 
end ; 

function phidot_o(t : float) s float; 
var tdo : float; 

U>gl n 

tdo := the t ado t_o ( t ) ; 
phidot_o : = tdo- (m*w_r ) ; 
end : 

function w_l ( t : float) : float: 

var betat» _ oo. pj. pdo> pdj, numer. denom : float; 
begin 

betat :* betaCt); 
oo : s phl_o ( t ) ; 
pj *. = phl_j(t): 
pdo := phTdot_o(t): 

Pdj : = phi do t_ j C t ) : 

numer be t a t * be ta t «pd j ♦odo* ( pd j ♦ odo ) *“be t a t *c os Cp j-po ) ; 

denom :* l*2'«betet*cos(oJ-po)*betBt , ‘betet; 
w_l :* m»*v_r*(numer/denom): 
end! 



<SF-> 



{ 5 F ♦ } 

function TNTargetFCt ; float) : float: 



{- This Is the function to Integrate -) 

{ ) 



var tempi : float; 
oe gl n 

tempi : = w_l ( t ) : 

TNTargetF := templ»r*cos(w_r*t ) ; 
end: { function TNTergetF ) 

{ 5F- ) 



{ CF • } 

function TNT ar ge tF_s 1 n ( t : float) ; float; 



{ > 

{- This Is the function to Integrate -) 

{ ) 



var tempi : float; 
begin 

tempi ;= w_l ( t ) : 

TNTarge tF_sln := tempi «r«sln(w_r»t ) ; 
end: ~ { function TNTarge tF_sl n ) 

(OF-) 



procedure 1 nl t 1 allze ( var LowerLlmit : float; 

var UpperLlmlt : float: 
var Integral ; float; 
var Tolerance : float; 
var Mavintervals : Integer; 
var Numlrtervals : Integer; 
var Error : byte); 



{ } 

{- Output: LowerLlmit. UoPterLlmlt. integral, Tolerance. -) 
{- Max Intervals . Numlntervais . Error -) 
{- -) 
{- This procedure Initializes the ar.ove variables to zero -) 
< ) 



begin 
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wrlteln; 

LowerLlmlt := 0: 

Upper-Limit !* 0; 

Integral :* 0; 

Tolerance ;* 0: 

Max Intervals := 0; 

Numlntervals := 0; 

Error := 0: 

end: { procedure Initialize ) 



procedure Re_I nl 1 1 a 11 ze ( 
var Integral : float: 
var Numlntervals : Integer; 
var Error : byte); 

begin 

Integral ; = 0; 

Numlntervals r= 0; 

Error ;= 0; 

end; { procedure Reinitialize > 

procedure get_root_name_of_f lies ( var Filename : string; 

def ault_f llename : string); 

begin 

wrlteln: 

FlleName := de f aul t_f 1 1 ename ; 

wrltet'Enter root of file name t’» FlleName. ’) '); 
r e adln (Filename ) : 

If Filename = * ' then FlleName := de f aul t_f 1 lename ; 
end; { procedure GetOutputFlle_Name > 



procedure GetOutputFl le_Name ( var OutFlle : text; 

var FlleName ; string; 
defaults llename : string); 

var 

Ch : char; 

begin 

repeat 

Ch := * V * ; 
wrlteln: 

FlleName := default_fllename: 
write(’Enter file name I*. FlleName* ’} ’); 
readln(FHeName ) : 

If FlleName = ' ' then FlleName :* de f aul t_ f 1 1 en& ne ; 
assl 9 n t Ou tFl 1 e » FlleName); 
reset! OutFlle); 

If loresult = 0 then { The file already exists. > 

begin 

cl ose ( OutFl 1 e ) ; 
wrlteln: 

write! 'This file already exists. ’); 
write! 'Write over It (Y/N)? ’); 

Ch := uocase (ReadKey ) ; 
wr lteln(Ch): 

If Ch = 'V' then assl gn (OutF) ] e , FlleName); 
end : 

If Ch = 'Y* then 
begin 

rewrl te (OutFlle ) ; 

IOCheck ; 
end ; 

until ((Ch * ’Y') and not(IOerr)); 
end: { procedure GetOutputFl le_Name ) 

procedure GetDatafvar LowerLlmlt : float; 

var Upper Llml t : float; 

var Tolerance : float; 

var Maxlntervals : Integer); 



{- - } 

{- Output: LowerLlmlt. UpperLlmlt. Tolerance. Maxlntervals -> 
(- -> 
{- This procedure assigns values to the above variables -) 
{*- from keyboard Input -> 
( - ) 



procedure Ge tLlml ts ( var LowerLlmlt : float: 
var UpperLlmlt : float): 



( > 

(- Output; LowerLlmlt. UpperLlmlt -> 

(- -> 

This procedure assigns values to the limits of -} 

(- Integration from keyboard Input -} 

( } 



begin 

repeat 

repeat 

LowerLlmlt := 0.0: 

IOCheck ; 

until not IOerr: 

wrlteln; 

repeat 

UpperLlmlt := 2«pi/w r; 

IOCheck; 

until not IOerr; 

If LowerLlmlt = UpperLlmlt then 
begin 

wrlteln; 

wrlteln(' The limits of Integration must be different.* 

wrlteln; 
end ; 

until LowerLlmlt <> UpperLlmlt: 
end; { procedure GetLl* Its) 



procedure GetToleranceCvar Tolerance : float); 

< } 

{- Output: Tolerance -> 
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{- This procedure reads In the accepted Tolerance -) 
{- from the keyboard. -> 



const def ault_tolerance = le~2; 
be 9 ln 

wr 1 te In ; 
repeat 

Tolerance s= default^tolerance ; 

wrl te ( 'Tolerance In answer: (> 0): ’): 

ReadFloet (Tolerance); 
lOCheck : 

If Tolerance <= 0 then 
begin 

IOerr := true; 

Tolerance := def ault_tolerance ; 
end ; 

until not IOerr; 

end; { procedure GetTolerance } 

procedure GetMaxIntervalstvar Maxlntervals : Integer); 



{- Output: Maxlntervals -) 
(- -> 
{- This procedure reads In the maximum number of -) 
{- subintervals to be used In approximating the -) 
{- Integral. Input Is from the keyboard. -) 
{ - ) 



be 9 l n 

wr i t e In : 
repeat 

Maxlntervals 1000; 

wr 1 1 e (’ Maxi mum number of subintervals (> 0): ’); 
Readlntt Maxlntervals); 

IOCheck : 

If Maxlntervals <= 0 then 
begin 

IOerr := true; 

Maxlntervals := 1000: 
end : 

until not IOerr: 

end; { procedure GetMaxI ntervals ) 

begin { procedure Ge Data ) 

GetLlmltsCLowerLlmlt. UpperLlmi t ) ; 

GetTolerance(Tolerance): 

GetMavIntervals(MaxIntervals); 

{ GetOutPutFile (OutFi le ) ; ) 
end; { procedure GetData ) 

procedure ResultstLowerLlml t : float: 

UpperLlmit : float: 

Tolerance : float; 

Maxlntervals : Integer; 

Integral : float; 

NumI ntervals : integer; 

Error : byte); 



{ > 

{- This procedure outputs the results to the device OutFlle -} 

{ > 

begin 

wr 1 1 e 1 n ( Out FI 1 e ) ; 

wrl teln(OutFile . ’Lower LimJt^SS* LowerLlml t :25 ) : 
wri teln (OutFi 1 e . ’Upper Llmlt:’:35. Uppe rLlmi t :25 ) ; 
writeln(OutFlle . ’ Tol er anc e : ’ : 35 . Tolerance :25 ) ; 

wr 1 1 eln (Out Fi le . ’Maximum number of subl ntervals :*: 35 » Max I nt ervals : 5 ) ; 
wr l t eln (Ou tFl 1 e . ’Number of subintervals used:’:35. Numlntervals :5 ) ; 
if Error = 3 then 
Displaywarning ; 

1 f Error In [ 1 . 2 ] then 



DlsplayError ; 



case Error of 

0 : wrl telnCOutFi le . * Integral :’ :25 . Integral); 

1 : wri teln(OutFl le • ’The tolerance must be greater than zero.*): 

2 : wr 1 1 e 1 n (OutFi l e . 

’The maximum number of Intervals must be greater than zero. 

3 : begin 

wr 1 te In (OutFi le . ’The Integral was not found with ’* Numlntervals 
’ subintervals . ’ ) ; 

wrl teln(OutFl le . ’The Integral thus far: ’, Integral): 
end : 

end; { case ) 

end; { procedure Results ) 

procedure 9 e t_be t a_ j ( var beta_J ; float); 
begin 

wrl teln; 
repeat 

beta_3 := beta_i_def ault ; 
wrl te ( ’ beta_ J : ’ ) ; 

ReadFloat (beta 1); 

IOCheck: 

If beta_3 < 0 then 
begin 

IOerr := true: 
beta_J := beta_J_de fault : 
end : 

until not IOerr; 

end; { procedure get_beta_j ) 

procedure 9 e t_del_theta (var del_theta : float); 
begin 

wrl teln ; 
reoeat 

del_theta := del_theta_de faul t : 
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write( ’del _t beta *• ' ) : 

Read c loat (del_theta ) ; 

IOCheck: 

1 f de l_the ta < 0 then 
be9i n 

IOerr :s true: 

del_theta := de l_thetB_def aul t : 
end : 

until not IOerr; 

end; { procedure 9 et_del_thete) 

procedure get_rho_3 ( var rho_J : float); 
begin 

wr 1 1 e In : 
repeat 

r h o_ 3 := rho_3_def suit ; 
wr i te ( ’ rho_3 : ' ) : 

Read r loat (rho 3 ) : 

IOCheck : 

If rho_3 < 0 then 
be 9 ln 

IOe r r : s true; 
rho_3 := r ho_3_ d e f au 1 t ; 
end ; 

until not IOerr: 
rho_3_start : = rho_3 : 

end; ~ { procedure ge t_de l_the t a > 

procedure ge t_the t bc_ 3 ( var thetac_3 *• floatJ; 
begin 

w r 1 1 e 1 n ; 
repeat 

tnetac_3 := t he t ac_ 3_d e f a ul t : 
wri-reC'thetac,! : *); 

Pe idFlcat ( thetac_3 ) ; 

ICIhec k ; 

1 5 the tac_3 < C then 
u e g 1 n 

IOerr : = true: 

tnetac_3 := thetac_3_default: 
end; 

until not IOe r r; 

end; { procedure 9 et_del_theta > 

procedure 9 et_del_w_r_rev (var de)_w_r_rgv : float): 
be 9 l n 

w r 1 1 e 1 r> ; 
repeat 

del e v •. = del_w_r_re v_def ault : 
write?’ de }_u_r_rev : 

PeadFloat ( del_w_r_r e v ) ; 

IOCheck; 

1* del_w_r_rev < 0 then 
begin ” 

IOerr := true: 

d ei_w_r_r e j := de l_w_r_r e v_de fault ; 
end : 

until not IOerr: 

eno; { procedure 9 e t_del_the ta > 



procedure lnlt_pa^ams; 
begin 

ge t_te t a_3 ( be t a_ j ) : 
ge del_?heta(del_ theta): 
ge t_r ho.. 3 ( r hc_ j ) 

9 e t_the tac_3 (thetac_3): 
get_del_w_r_rev Cael_w_r_rey ) ; 

w_r := w_r _r e u » 0 “ P 1 : 

del _w_r :=~eel_w_r_r ev“C*pl : 

t_c y c : - Z-Pi/*_r; 
w_r _ 3 :« (w.r*Cel>.r); 

t_C'/C_3 -• = Z-Cl/V. r •) ; 
ref _t := thet£r_.j/w -_j; 
oel_t_Z := del. the t e '2/ w_r_j : 
u_I :m := re*_t*oel_t.Z; 

while u_lir> • t_c do u_lim : = u_lim- t_cyc_3 ; 
l_i i fi r - r e ‘_t -de l_t_Z : 

w lie 1 Ilf*' < 0 do~i_i:m : = 1_1 1 rr,* t_cyc_ j : 

rho : = rh o_i r,3 t : 
thetac_T := thetac_T_J ni t ; 

end : 



procedure wr l t e_on e_l t e r a t 1 on_l 1 ne ( v ar out_file : text): 
begin 

wr 1 te In ( out_f 1 1 e » 

’ ’ . rho 3:6:3. 

* ’ . r h o ? 7 : 2 • 

* thetac_T:7:3 . 

*>'. x tj:~:2. 

‘ * . vlt3 :?:S. 

’ ’ • ran S e_*gt :?:3 . 

* ' • theta t s t : 7 : 3 . 

* • , reuj :7:3- 
’ * . n e w_y : 7 : 3 ) ; 

end ; 



procedure cs3_seeker ( time : float; 

var seeker_theta_prev . 
t ?eke"_pr.l o r ev. 
seer.er_theta._cur. 
seer, e r _ph j_* ur ; float); 
var s p , y p . v c . yc . x d . yd T float: 
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Begin 

UpperLlmit : = time: 

LowerLlmlt :* Upp e rLl m 1 t- t.c yc : 

> p : = s e e x e r _t h e t a_pr e v - f oc s 1 .1 eng t h ; 
yp :• seeker~ph'_prev»f ocal_i?ngth : 
vc ;= seekerltheta.cur * focal.length ; 
yc : = seeker~phl_cur>*f ocal. length; 
yd : = yp-yc:“ 
yd : = yo-yc: 

rho := sart'(yd“yd)*(yd*yd)); 
thetac.T : s atanC'yd. yd); 

» dr Dtl ve.Gauss.Quadrature(LowerLlml t . Upper-Limit. Tolerance. May I nt ervals . 

Integral. Num 1 n t er va 1 s . Error, STNTargetr); 

x. t j . I n t egral/ ( pi *m ) : 

Pe.l r>l 1 1 al lze C I nt e gra 1 . NumI nt ervals . Error); 

Adaptl ve.Gauss.Quaar atur e (LowerLlmi t . UpperLlmit, Tolerance, Maxlntervals . 

Integral, NumI nt ervals . Error, STNTar 9 e tF_sl n ) ; 

y. tj := lntegral/(pl*m); 

{ renge_t 9 t:=$ort(y_tj“<_t 3 *y_t J "y.t 3 ) ; 
theta_Tgt;=atan2(y”tl,y_t3): 

Oel.range:=range.tgt: 
del~x:=oel_ranoe“ccs(theta_tgt): 
del_y:=dellrange“sln( theta_tgt ) ; ) 

del_y : = x_t j ; 
del_y ; = y_t3; 

yc :a vp-del.x: 
yc := yp-cel.v: 

seeke»"_tneta_c ur -.= yc/ f oc el.length: 
seey er“pni_cur := yc / f oc al.i eng t h ; 
seeker_t M e • a_rrev : = seeke?_thets_eur ; 
seeker”phl_pr ev ;* seeker_pni_cur ; 
end : 

procedure wr 1 t e_t afcle.hpader ( var out. file : text); 
begin 

wr 1 1 e 1 n ( out_F 1 1 e ) ; 
wr l t eln( out"f 1 le . 

* 1 • . 

’ r ho_3 ’ . 

’ r ho * . 

' tht T* , 

* y tr . 

' y_tr. 

* r n g_t g ' , 

* tM~tg' . 

’ np w_x * . 



procedure wr 1 1 e.end.l t er e tl ons_l lne ( van out_file : text); 
begin 

wr l tel n ( ou t.f 1 1 e . 

1:6. 

’ ’ • r h 0 _ j : 6 : 3 . 

' • , rho : 7 : 3 , 

* tnetac_T;7;3, 

*> End . . . ’ )T 

end ; 

procedure wr i t e.f 1 1 e.he ader : 
begin 



wrlteln(csj 


Ml?): 








wrltelMcsj 


file. 


* rr = ’ . m : C ) 






wr: * ein( cs j 


file. 


’ = ’ » m_3 : 


6 ); 




wr 1 1 eln t c s : 


file. 


* r = • , r : 7 : 


3 ) ; 




writelnfcr) 


f lie . 


* beta. 3 


= • . 


beta.3 : 7 : 3 ) ; 


w’-lteln(esj 


“file. 


• w r rev 


= ’ , 


w.r.rev :C ) : 


wrltfiln(cs: 


“file. 


•del w r rev 


- ’ , 


dpl~w_r rev : 7: 3 ) : 


wrJtelntCSl 


file. 


’ p n 1 c_r 


= 1 , 


p h 1 c _r : 7 : 3 ) ; 


w r 1 ♦. e 1 n C c S 3 


.file. 


’ Phi 


= ’ , 


phjc_3:7: 3): 


wrl-elncsj 


e i j e . 


' rho_Tr l t 


= ■ , 


r ho_Inl t : 7 : 3 ) ; 


writeln(csj 


file. 


’ rbO_3_c tart 


= ’ . 


rho_j_st?r t ; 7 : 3 ) ; 


w-~ 1 1 e 1 n ( c S 3 


Mie. 


’d-'l.theta 


s 1 , 


del.t het a : 7 : 5 ) : 


wrltelntcsj 


If J le . 


* t he T ac_T_l nit 


s • , 


thetac_T_lni t : 7 ; 3 ) 


writelntcsj 


file. 


’ ♦ he * a 3 


= * . 


theta. 3:7:' > ; 


write 1 n ( c S 3 


“file. 


• f 1 1 :_rO"S 


= • , 


flit_ cons: 7:3): 


wrltelntcsl 


"file. 


’ f ocal.iength 


= * . 


focal. length: 7:3); 


wrltelnCcs? 


f ; 1 p ) ; 








wrlteln(cs3 


f lie ) ; 








wr 1 1 e_t ao 1 e 


he?oer(csi file); 






wr i telnl c sj 


_f 1 1 e ); 









end ; 



begin { main program ) 

Cl rC-c r ; 

get_root_nane_r'^_flles(dp<ault_< : llen3me. 1 nl t.de f aul t_ filename); 

Ge tOu tD’.-TF 1 le~N»me ( c s J.f 1 le • FlleName. default_fllename 4 ’cs3’ ): 

1 nl t .oar amo : 

Ini tlallze ( Lc..* rLlml t . UpperLlmit. Integral. Tolerance. Maxlntervals. 
NumI nterva) S • Error): 

GetData (LowerLlni t • UrperLiMt. Tolerance. Maxlntervals); 

range_tgt ;= 0; 
the * a.t gt ;t o : 
r •’P : = r nc.i nit; 
tnetac_T t he tac.T.l ni t : 
w- l te_~l le.he ader 



end . 



{ unit cslunit) 
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